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Abstract 

We calculate the leading-order matrix element for the decay B 
K*"f in the quenched approximation of lattice QCD on a 24"^ x 48 lattice 
at /3 = 6.2, using an 0(a)-improved fermion action. Extrapolating the 
quark masses to their physical values we obtain an on-shell form factor of 
ri(g2=0) = O.I5II4, where the errors quoted are purely statistical. We 
find Ti is approximately independent of the spectator quark mass and 
extract Ti{q^=0) = 0.151*^4 if this independence is assumed. We compare 
this with the same form factor derived (in the Standard Model) from the 
CLEO experimental branching ratio of BR{B K*^) = (4.5 ± 1.5 ± 
0.9) X 10~^ and find the results to be consistent within statistical errors. 

To be submitted to Phys. Rev. Lett. 



PACS: 12. 38. Go, 13.40.Hq, 14.40.Jz 



^ inf o-ukqcdOhep.ph. soton. ac . uk 



B K*^ : Penguins on the Lattice 

UKQCD Collaboration 
February 1, 2008 

Abstract 

We calculate the leading-order matrix element for the decay B 
K'j in the quenched approximation of lattice QCD on a 24^ x 48 lattice 
at (3 = 6.2, using an 0(a)-improved fermion action. Extrapolating the 
quark masses to their physical values we obtain an on-shell form factor of 
ri(g2=0) = 0.15114, where the errors quoted are purely statistical. We 
find Ti is approximately independent of the spectator quark mass and 
extract Ti{q^=0) = 0.15^4 if this independence is assumed. We compare 
this with the same form factor derived (in the Standard Model) from the 
CLEO experimental branching ratio of BR{B K*'y) = (4.5 ± 1.5 ± 
0.9) X 10"^ and find the results to be consistent within statistical errors. 

1 Introduction 

Theoretical interest in the rare decay B — > K*^ as a test of the Standard 
Model has recently been renewed by the experimental results of the CLEO 
collaboration For the first time, this mode has been positively identified 
and a preliminary determination of its branching ratio given. 

The significance of i? — > K*^ arises from the underlying flavor-changing 
quark-level process h — > 57, which first occurs through penguin-type diagrams 
at one-loop in the Standard Model. The process is also sensitive to new physics 
appearing as virtual particles in the internal loops. Existing bounds on the 
6 — > 57 branching ratio have been used to place constraints on supersymmetry 
[H H || and other extensions of the Standard Model ||, H . 

In order to compare the experimental branching ratio with a theoretical pre- 
diction it is necessary to know the relevant hadronic matrix elements. These 
have been estimated using a wide range of methods, including relativistic and 
nonrelativistic quark models 0, |[ |^ , two-point and three-point QCD sum rules 
[ p^ [TTl [1^ , [1^ and heavy quark symmetry |Q, but there remains some dis- 
agreement between the different results. It is therefore of interest to perform 
a direct calculation of the matrix elements using lattice QCD. The viability of 
the lattice approach was first demonstrated by the work of Bernard, Hsieh and 
Soni H in 1991. 
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In the leading-log approximation the B — )■ K*j transition is caused by a 
single chiral magnetic moment operator from the effective weak Hamiltonian. 
In the notation of Grinstein, Springer and Wise this is, 

07 = Y|^mfcsa^,i(l + 75)5i^^^ (1) 

with an on-shell matrix element given by, 

M = ^^0C^{m,)Vt,V;,v'*{K*\sa,,q-^bn\B), (2) 

where q and rj are the momentum and polarization of the emitted photon. The 
coefficient CT{mb) arises from the mixing of O7 with other effective operators 
in running the scale down from Mw to mf,. The anomalous dimension matrix 
of all the effective operators at the one-loop level has been calculated by several 
groups and is now well- understood fl^ . 

Following Bernard et al. the general matrix element can be parametrized in 
terms of the momentum, k, and polarization, e, of the K* ^ and the momentum, 
p, of the B meson, using three form factors, Ti, T2 and T3, where Ti is chosen 
to be real, so that T2 and T3 are purely imaginary, 

{K*{k,e)\J,\B{p))=ClT,{q^) + ClT2{q^)+ClT:,{q^) (3) 

Jfi ^su^^q^bn, q = p-k, (4) 

Cl = 2e^,ApeVfc'' (5) 
= ef,{m%-mj^,)~e-q{p + k)f, (6) 

Cl = e.q{q,--^l-^ip + k),). (7) 

The on-shell (g^ = 0) matrix element depends on Ti only, since T2{q^=0) = 
— «Ti(g^=0) and the coefficient of T3 is zero. Performing the necessary phase 
space integral and sums over polarization vectors gives the decay width for 
B^K*j, 

r{B ^ K*^) = ^mlGlml (l - ^) \Vt,V;f\Cr{m,)\'\T,{q^^O)\\ 

(8) 

By computing the matrix elements on the lattice for various q^, the on-shell 
value of the form factor Ti (0) can be obtained by interpolation. 



2 Computational Details 

We work in the quenched approximation on a 24^ x 48 lattice a.t /3 — 6.2, 
which corresponds to an inverse lattice spacing a^^ = 2.73(5) GeV, evaluated 



2 



by measuring the string tension Pq |. Our calculation is performed on sixty 
5/7(3) gauge field configurations (for details see Refs. Q and Q). The quark 
propagators are calculated using an 0(a)-improved Wilson fcrmion action ]l9| . 
We use gauge-invariant smeared sources for the heavy quark propagators with 
an r.m.s. smearing radius of 5.2 [Q. Local sources are used for the light quark 
propagators. 

As the mass of the b quark is almost twice the inverse lattice spacing, direct 
computation of a 6-quark propagator is not feasible. We therefore compute 
heavy-quark propagators with masses in the region of the charm-quark mass, 
and extrapolate. 

Our statistical errors are calculated according to the bootstrap procedure 
described in Ref. iQ, using 250 bootstrap samples. 

To obtain the matrix element (y(fc)|so'p^6|P(p)), we calculate a ratio of 
three-point and two-point correlators, 

„ ,, , . ^ CXJf.tf'i'-il im 

'^""■"■'''^ = c-(„-,,p)c;"(,p-a' 



where. 



x,y 



t,tf—t~too^ — ' 2Ep 2Ev 



and 



Cr{t,p)= E^e'^^-(4(i,f)Jp(0)) ^ ^e""""* 

CT{t,k)^ -(0E5e*-74„(t,f)J^(O)) ^-.^^e-^-* (11) 

with Jp and Jy interpolating fields for the pseudoscalar and vector mesons re- 
spectively. T^y is the 0(a)-improved version of the operator baf^^s [pl| . The full 
matrix elements can then be derived by using the relation a^^j^ — —■^^^uXpCr^'^- 
We employ time reversal symmetry to obtain the correctly ordered matrix el- 
ement, ly {k)\'saaiMP{p)) . To evaluate these correlators, we use the standard 
source method |2^. We choose t/ = 24 and symmetrize the correlators about 
that point using Euclidean time reversal [^. We evaluate Cp^j^ for three 
values of the light quark mass {m = 0.14144,0.14226,0.14262), two values of 
the strange quark mass {k,s = 0.14144,0.14226) which straddle the physical 
value (given by nP^y = 0.1419(1) (ID), and two values of the heavy quark 
mass (/t/i = 0.121,0.129). We employ two values of the B meson momentum 
((12a/7r)p = (0,0,0), (1,0,0)), and seventeen values of the momentum, g, in- 
jected at the operator, with magnitudes between and 27r/(12a). To improve 
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statistics we average over all equivalent momenta, and utilise the discrete sym- 
metries C and P, where possible. 

Provided the three points in the correlators of Eq. ( ^) are sufhciently sepa- 
rated in time, the ground state contribution to the ratio dominates: 

^"'^'^ M.^oo E e)|s'T^^6|F(p)) + . . . (12) 

and Cp^i, approaches a plateau. The factors Zp, Zy and the energies of the 
pseudoscalar and vector particles are obtained from fits to two-point Euclidean 
correlators. 

The form factor Ti can be conveniently extracted from the matrix elements 
by considering different components of the relation, 

A{k"pP - p'^k^)T,{q') ^ e"'''"^Y.'p{V{k,e)\sa^,b\P{p))q''. (13) 

e 

We see a plateau in Ti about t — 12, and fit Ti{t;p,q) to a constant for t = 
11, 12, 13, where correlations are maintained between all of the time slices. The 
use of smeared operators for the heavy quarks provides a very clean signal, with 
stable plateaus forming before timeslice 11. Data with initial or final momenta 
greater than (7r/12a)V2 are excluded, because they have larger statistical and 
systematic uncertainties. 

The data for the heaviest of our light quarks, ki — Kg — 0.14144, with the 
smallest statistical errors, are shown in Fig. nl 



3 Results 

We fit Ti{q'^) to a linear model in order to obtain the on-shcU form factor 
Ti(g2=0), 

T^{q^) = a + bq^, (14) 

allowing for correlations between the energies of the vector and pseudoscalar 
particles and Ti at each q^ . For our range of masses and momenta the differences 
between linear and pole model fits are small. The data and fit for ki — Kg = 
0.14144 are shown in Fig. |[ 

The light quark mass is set to zero by a correlated chiral extrapolation to 
Ki = Kcrit(= 0.14315(2)) We assume that the on-shell Ti varies with the light 
kappa values, ki, according to a linear model, 

Tl{Ks,Kh, m) = T^"\Ks,Kh) + Ai{Ks, Kh) (- —j ■ (15) 

The strange quark mass is set to its physical value by interpolation (k^ = 
i^phys ^ 0.14183(5)). 
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Finally we perform an extrapolation from the two pseudoscalar masses up 
to ms using a model motivated by heavy quark effective theory, 

Ti(g2^0;mp) = A+— . (16) 

■nip 

After performing this extrapolation, we obtain Ti{q^^{)]mB) = 0.15^]^4, 
where the quoted error is purely statistical. The finite renormalization needed 
for the lattice-continuum matching of the a^j^y operator has been calculated [ p5[ 
but has a neghgible effect here (0(2%)) and is not included. 

We note that the slopes of the form factor Ti with respect to ni in the 
chiral extrapolations are consistent with zero (Fig. which indicates that 
Ti{ks, Kh,Ki) is almost independent of k;. However, this behaviour occurs only 
for the spectator quark, and is not seen to the same extent for the interacting 
strange quark. We explore this by fitting Ti to a constant for the three values of 
Ki. We find that the per degree of freedom is comparable to the original linear 
model, indicating that the model is statistically valid. Using this approach, the 
final statistical error is significantly reduced, and we obtain Ti(g^=0;TOB) = 
0.15i^. The results for Ti, using both analysis procedures, are shown in Fig. W. 

4 Conclusions 

In this letter we have reported on an ah initio computation of the form factor 
for the decay B K*^. The large number of gauge configurations used in this 
calculation enables an extrapolation to the appropriate masses to be made and 
gives a statistically meaningful result. To compare this result with experiment 
we convert the preliminary branching ratio from CLEO, BR{B K*"f) = 
(4.5 ± 1.5 ± 0.9) X 10~^ based on 13 events, into its corresponding Ti form 
factor, assuming the Standard Model. We work at the scale fi = mt, = 4.39 GeV 
and use values from the Particle Data Book ||2^ combined with Eq.( ||). Setting 
the mass of the top quark to be rrit = 100, 150 and 200 GeV we find T'^^^ to be 
0.23(6), 0.21(5) and 0.19(5) respectively. The two lattice results are consistent 
with these experimental numbers within statistical errors. This is also shown in 
Fig. |. 

Although systematic errors of this calculation resulting from the quenched 
approximation, finite volume and other lattice artefacts remain to be explored, 
we believe that we have shown the phenomenological utility of the lattice for 
probing the limits of the Standard Model. 
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Figure 1: Ti vs. timeslice, t. (For computational reasons, only timeslices 7 — 16 
were stored) 
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Figure 2: Ti{q^), with a linear fit. The dotted Hnes represent the 68% confidence 
levels of the fit at = 
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Figure 3: Chiral extrapolation of Ti(g^=0). The dotted lines indicate the 68% 
confidence levels of the fit. m, Hght is the lattice pole mass. 
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Figure 4: Extrapolation of Ti(g^=0) to tob- LCE — using linear chiral extrap- 
olation, CCE — using constant chiral extrapolation for the spectator quark. 
(N.B. for clarity, the LCE and CCE points have been displaced horizontally by 
0.02 to the left and right respectively) 
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